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SIGNIFICANCE  AND  EXPLANATION 

In  experimental  work  it  is  sometimes  desired  to  make  n 
different  inferences  from  a given  data  set.  Associated  with 
the  data  is  a hypothesis  called  the  overall  null  hypothesis 
which  is  to  be  tested,  and  under  the  assumption  that  the  overall 
null  hypothesis  is  true  it  is  desired  to  have  a preassigned  high 
probability  1-a  of  having  all  the  individual  inferences 
be  correct.  If  the  proper  assumptions  are  made  concerning  the 
statistical  distributions  involved,  this  requires  finding  the 
abscissa  t^2n  ^or  a sPecifie(*  t-distribution  having  area 
ct/(2n)  to  the  right  of  it.  Such  abscissas  are  called  Bonferroni 
t-percentiles . Since  statistical  tables  only  give  abscissas 
corresponding  to  a,  it  is  sometimes  troublesome  to  find 
ta/2n*  Whereas  previous  work  of  the  author  is  applicable  if 
a/(2n)  is  small,  the  present  paper  gives  a method  that  works 
well  for  any  a/(2n)  between  .0  and  .5  and  is  readily 
implemented  on  a computer. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


Rational  Bounds  for  the  t-Tail  Area 


With  an  Application  to  Bonferroni  t-Percentiles 
Andrew  P.  Soms 

The  purpose  of  this  paper  is  to  complete  the  extension  to  the  t- 
distribution  of  results  known  for  the  normal  (see,  e.g. , Johnson  and  Kotz 
1970,  Chapter  33).  In  Soms  (1976)  an  asymptotic  expansion  for  the  tail  area 
of  the  t-distribution  was  given  good  for  small  areas  and  in  Soms  (1977)  some 
improved  bounds  were  given,  none  of  which,  however,  had  simultaneously  the 
right  Uniting  behavior  as  x -*•  0 or  x ®.  In  this  paper  lower  and  upper 
bounds  are  given  with  this  property,  the  accuracy  of  the  approximation  is 
studied,  and  the  results  applied  to  the  calculation  of  Bonferroni  t-percen- 
tile  and  descriptive  levels. 

2.  Extensions  of  Boyd's  Results 

—1/2 

Boyd  (1959)  showed  that  if  4>  (x)  = (2x)  e-X  ' 2 , F^(x)  ■ /°°  $(t)dt, 

x 

and  R^  * F^(x)/4i(x) , x > 0 , then 

p<x’Y»io>  * R* 

2 2 2 1 /2 

where  p(x,y)  = (y+l)/[(x  + -(y+1)  ) ' + yxl , y ■ 2/(ir-2),  y , - ir-1, 

n max  min 

and  the  bounds  are  the  best  possible  in  the  class  {p(x,y),  y > -1).  This  is 
also  discussed  in  Johnson  and  Kotz  (1970,  Chapter  33).  Here  we  will  obtain 
analogous  results  for  the  t-distribution. 


Let 


, -(k+l)/2 

fk(t)  = ck(l  + tVk) 


k an  integer  >1,  and 


B r((ktp/2) 
r(k/2)(irk)Js 


Fk(x)  = l~Fk(x) 


fk(t)dt 


9 


for  x > 0 her  throughout  the  paper.  Also  let 
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R = F (x)/((l+x  /k)f.(x) ) 

x k k 


was  shown  in  Soms  (1976)  that 


1/x  - k/((k+2)x  ) < Rx  < 1/x 


Rx  = 1/x  - k/((k+2)x3)  + 0(x'5) 


(2.1) 


In  Soms  (1977)  it  was  shown  that 


(k-l)x 


+ 1 + 


(k+l)x}: 

2k  J 


< R < 

x 4k 


k-1  [’(k+Dx)2)**" 

2k  { 4k  J J ’ 


(2.2) 


where  the  lower  bound  is  valid  for  k > 1 and  the  upper  for  k > 2.  It  can 
be  seen  that  these  bounds  do  not  have  the  right  limiting  behavior  as  x + 0, 
Motivated  by  (2.2),  consider  the  problem  of  finding  best  upper  and  lower 
bounds  of  the  form 


(x  +BP  + yx 

Since  lim  R = l/2c,  and  R ~ 1/x  (as  x -*•  °°),  one  must  have 
x k x 

x-*-0 

a = 1 + y and  a/B  = i/2ck  • 

2 2 

Thus  B = 4c^  (1+y)  and  the  approximating  functions  p(x,y)  must  be  of  the 


/ x 1+Y 

p(x,Y)  = — 5 5 o-r 

( x +4ck  Z ( 1+y ) ^ )^+Yx 

For  small  x, 


R = l/2c,  - x + 0(x  ) 
x k 


-2- 


and 


p(x,y)  = l/2ck  - Yx/(4ck2(l*Y))  ♦ 0(x2)  . 

For  large  x,  from  (2.1), 


Rx  = 1/x  - k/((k+2)xJ)  ♦ 0(x"5) 


and 


p(x,Y)  = 1/x  - 2c /(ltY)/x3  ♦ 0(x'5)  . 


So  if  p(x,Y ) is  to  maximize  R^  it  is  necessary  that 

Y 


4ck  (y+1) 


-1  and  > 2c  2(1+y)  , 


or 


4c, 


-1  < Y < 


l-4c. 


and  -1  < y < 


-1 


2(k+2)c, 


or 


-1  < Y < min 


4C, 


[l-4c  2(k+2)c, 


- 1 


(2.3) 


Note  that  from  Lemmas  A2  and  A3  and  the  direct  evaluation  of  — 


- 1 


2(k+2)c. 


for  k = 1 the  interval  in  (2.3)  is  non-empty.  Similarly,  if  p(x,Y)  is  to 


minimize  R^,  it  is  necessary  that 


4ck  (y+1) 


>1  and  jJj  < 2ck2(Y*l)  , 


or 


l-4ck2  < 2(k+2)ck2 


is  equivalent  to 

ck2  “ k/(6k+4)  . 

Hence  from  Lemma  A4,  if  k > 3,  then  for  maximizing  it  is  necessary  that 

Y 1 4ck2/(l-4ck2)  , 
and  for  minimizing 

Y > j - 1 * 

2(k+2)ck 

If  k = 2,  then  4ck2/(l-4ck2)  - k/(2(k+2)ck2)  - 1,  and  it  is  verified  by 
differentiation  that  (here  and  throughout  will  denote  differentiation) 

- |(1  + ^r-)fk(x)p(x,Y)j  = fk(x)  , 

2 2 2 

and  hence  in  this  case  = p(x,y)»  where  Y = 4ck  /(l-4ck  ) = k/(2(k+2)ck  ) - 1. 
Therefore  for  k = 2 the  approximating  problem  has  been  solved  in  the  best 
possible  way  and  this  case  will  not  be  considered  further.  If  k = 1,  then 


Using  the  above  as  motivation,  the  first  result  is 
2 2 

Theorem  2.1:  Let  Y = 4c^  /(l-4ck  ).  Then  for  k > 3,  R^  < p(x,y)  and  for 
k = 1,  Rx  > p(x,y ) • 

2 

Proof:  Note  that  R'  = (.— R - l)/(l+x2/k)  and  R"  = ((JSll)(lk-l>.x.  +1)R 
X K X k k x 

(k  32x)/(1.t.x2/k).  Let  f ( x)  = p((x,Y)  - R • Then  f(0)  = f'(0)  = 0 and 
K x 

f"(0)  = ~±+2L — 5-  - 

8ck3(l+Y)2  2kCk 

16kck**  + (4k+4)ck^  - k 
8ck3(l+Y)2  k 

It  is  verified  directly  that  f"(0)  < 0 for  k = 1.  For  k > 3,  f"(0)  > 0 is 


equivalent  to 


16kck4  + (4k+4)ck^  > k , 


(2.4) 


which  will  certainly  be  true  if 


4 2 

16ck  + 4ck  > 1 . 


(2.5) 


In  Lemma  A1  it  is  shown  that  the  even  and  odd  ck's  for^  an  increasing 
sequence.  It  follows  by  direct  verification  that  (2.4)  is  true  for 
3 < k < 16  and  that  (2.5)  is  true  for  k = 17  and  k = 18  and  hence  (2.4)  is 
true  for  k > 3.  So  for  k > 3,  f"(0)  > 0.  Consider 

h(x)  = xf/(l+x2/k)  - f'  . 


By  some  algebra,  it  follows  that 


2 . T 2 2 2^  ~t  ^ o 9 o *5 

(1+x  /k)|^(x+4ck  (1+Y)  ) +YxJ  (x  +4Cj(  (1+Y)  ) U-4^2)  h(x) 


= xf 


■Oc2+4c  ,2/ ( 1-4C.  2>2)*+J  f 1 + 4 k 


,*k2<8ck2-1>  71— 7^77 

Mx  +4c,  / ( l-4c. 
k k 


Ck2  + 16ck4  I*  (2-6) 


-5- 


Suppose  now  that  f(x)  <_  0 for  some  x < xQ.  Then  since  f(x)  is  increasing  at 
the  origin,  it  must  have  a maximum  in  (0,xQ)  for  which  f > 0 and  f'  = 0, 
contradicting  (2.10).  If  f(x)  < 0 for  some  x > xQ,  then  f must  have  a 
minimum  (lim  f(x)  = 0)  for  which  f < 0 and  f'  =0,  contradicting  (2.11). 

X-HJO 

So  R < p(x,Y)  for  k > 3. 
x 

Consider  now  k = 1.  By  (2.7),  (2.9),  and  the  fact  that  f"(0)  < 0 for 

k = 1, 

h(x)  > 0 , x < xQ  , (2.12) 

h(x)  = 0 , x = xQ  , 


-6- 


and 


h(x)  < 0 , x > x0  . 


(2.13) 


Suppose  that  f(x)  > 0 for  some  x < xQ.  Then,  since  f is  decreasing  at  the 
origin  there  must  be  an  x < xQ  for  which  f < 0 and  f = 0,  contradicting 
(2.12).  If  f(x)  > 0 for  x > xQ,  there  must  be  an  x > xQ  for  which  f > 0 
and  f'  = 0,  contradicting  (2.13).  Hence  for  k * 1,  R > p(x,Y),  completing 
the  proof. 

Consider  now  the  problem  of  minimizing  R^  for  k > 3 and  maximizing 
for  k = 1.  Using  similar  reasoning  to  that  used  to  arrive  at  the  statement 
of  Theorem  2.1,  we  have 


Theorem  2.2;  Let  Y = 


2(k+2)c, 


- 1.  Then  for  k > 3,  R > p(x,Y)  and  for 


k = 1,  R < p(x,Y). 


Proof:  Let  g(x)  = Rx  - p(x,y).  Then  g(0)  = 0 and 


g' (0)  = - 1 + f j - l)/|>ck2  J 

'2(k+2)c,_  VI  K 2(k+2)ck^ 


(_6k-4)ck  +k 

= i 2 ’ 

4ck 

and  hence  by  Lemma  A5,  g'(0)  > 0 for  k > 3,  and  g'(0)  < 0 for  k = 1.  Consider 


h(x)  = xg/(l+x2/k)  - g’ 


After  some  algebra. 


= [-< 


h=  |-(x2+  4ck2(l+Y)2)  ((Y-l)x2 + (ltY)(-4ck2(l+Y)  ♦ Y)) 


<{( 


+ x^(Y-l)x2+(l+Y)(4ck2(l+Y)((l+l/k)Y-l+l/k)  - 1) 


}] 


/[((x2  + 4ck2(l+Y)2)^+  Yx)  (x‘  + 4ck2(l+Y)2)^(l+x2/k)J 


* 


or,  upon  some  further  algebraic  simplification. 


jg  2 lg  Ho 

((x2+4ck2(l+lr)2)  +Yx)  (x2+4ck2(l+Y)‘ ) (ltx2/k)((x2+4ck2(l+Y)‘ ) +x)h(x)/(l*Y)’ 


= -4ck2[x2(Y-l)  + (l+Y)(-4ck2(l+Y)  +Y) 


>] 


♦ x2|4ck2((l+l/k)Y+l/k)  - ll  + (x2+4ck2(ltY)2)x 


• 4ck  (<l+l/k)Y+l/k) 


-4- 


(2.1 4) 


Note  that 


4ck2((l+l/k)Y+l/k)  - 1 = 4ck2(y-l)  - 4ck2<U+l/k)Y+I/k)  ♦ 1 


(2.15) 

Hence,  using  (2.15),  the  right-hand  side  of  (2.14),  apart  from  the  constant 
term,  can  be  written  as 


2 2 2 2 H 

-cx  + cx(x  +4ck  (1+Y)  ) 


2 , 2.,  .2 
c 4ck  (1+Y) 


c[i+(i+4ck2(in)2/x2)i5]  ’ 


(2.16) 


2 2 •> 

where  c = 4ck  ((l+l/k)Y+l/k)  - 1.  Also,  -4ck  (1+Y)  + Y < 0 is  equivalent  to 


ck  > k/(6k+4),  and  hence  by  Lemma  A5,  -4ck  (1+Y)  + Y > 0 for  k > 3 and 


-4ck2(l+Y)  + Y < 0 for  k = 1.  Further,  4ck2( (l+l/k)Y+l/k)  - 1 > 0 is  equiva- 


lent to 


ck  < k/(4k+8)  , 


and  since 


k/(4k+8)  > k/(6k+4)  , k > 3 , 


by  Lemma  A5 


4c  ((l+l/k)Y+l/k)  - 1 > 0 , k > 3. 


It  is  verified  by  direct  evaluation  that  4ck  ( (l+l/k)Y+l/k)  -1  < 0 for  k = 1. 


1 


-R- 


2 2 

Consider  the  limit  2cck  (1+y)  of  (2.16)  as  x -*■  00 . For  k > 3, 


2cck2(l+y)2  > 4ck2(l+Y)(-4ck2(ln)+Y) 


(2.17) 


is  equivalent  to 


„ 2 (k ) (k+4 ) 

Ck  8(k+l)(k+2! 


(2.18) 


(2.18)  will  certainly  be  true  if 


2 _ k+4 

Ck  8(k+2)  ’ 


(2.19) 


and  it  is  verified  directly  that  (2.19)  holds  for  k = 9 and  k = 10  and  hence 
(2.17)  is  true  for  k > 9 by  Lemma  Al.  It  is  also  verified  directly  that 
(2.17)  holds  for  3 < k < 8 and  that  the  inequality  is  reversed  for  k = 1. 


Hence  it  follows  that  for  k > 3 


h(x)  < 0 , x < xQ  , 


(2.20) 


and  for  k = 1 


h(x)  = 0 , x = Xq  , 

h(x)  > 0 , x > xQ  , 


(2.21) 


h(x)  > 0 , x < xQ  , 


(2.22) 


h(x)  = 0 , x = xQ  , 


h(x)  < 0 , x > x. 


(2.23) 


Consider  k > 3 first  and  suppose  that  g(x)  < 0 for  some  x < xQ.  Then,  since 
g is  increasing  at  the  origin,  there  must  be  an  x in  (0,xQ)  such  that  g > 0 
and  g'  = 0,  which  contradicts  (2.20).  If  g(x)  < 0 for  some  x > xQ  then 
there  must  be  an  x > xQ  for  which  g < 0 and  g'  = 0 and  this  contradicts 

(2.21).  Hence  > p(x,y)  for  k 3.  If  k = 1 and  g(x)  > 0 for 
some  x < x , then,  since  g is  decreasing  at  the  origin,  there  must  be  an  x 


in  (0,xQ)  such  that  g(x)  < 0 and  g'(x)  = 0,  contradicting  (2.22).  If 


g(x)  > 0 for  some  x > xQ,  then  there  must  be  an  x > xQ  for  which  g > 0 and 

g'  =0,  contradicting  (2.23)  and  completing  the  proof  of  Theorem  2.2. 

We  summarize  the  results  of  Theorem  2.1  and  2.2  below  in  an  algebrai- 

2 2 

cally  equivalent  form.  For  k > 3,  let  y - 4c.  /(l-4c,  ) and 

" ma  x K k 


Y„,. 


2(k+2)c, 


- 1.  Then 


( ( x/  ( 1+Y„,.  „ ) ) 2+4c^  ^ 2 ) + Ym, „*/  ( 1+Y, 


min 


mn 


i„>]  X<  V 


+ y x/(1+y  ) 

max  max _ 


-1 


(2.24) 


If  k = 2,  then  y = y . = y and 

max  min 


Rx  = jj((x/(l+Y))2  + 4ck2)  +Yx/(l+Y)J 


For  k = 1,  (2.24)  holds  with  the  definitions  of  y and  y . interchanged. 

max  mm  & 

Since  9p(x,Y)/3Y  < 0 for  x > 0,  it  follows  that  for  k > 3,  y is  the 

max 

best  value  of  y for  which  R < p(x,Y)  for  all  x and  y . is  the  best 

x min 

value  for  which  Rx  > p(x,Y)  for  all  x.  A similar  statement  holds  for  k = 1 

with  the  definitions  of  y and  y . interchanged.  Equivalently,  the 

max  mm 

bounds  given  in  Theorems  2.1  and  2.2  are  uniformly  best  in  the  class 
considered. 

The  accuracy  of  the  lower  and  upper  bounds  are  related  to  the  maximum 

value  of  the  ratio  p(x,Y  )/p(x,Y  • )•  Call  the  ratio  r(x)  and  consider  the 

max  min 

maximum  value  of  r(x),  which  must  be  attained  for  some  positive  x at  which 
r'(x)  = 0.  By  differentiation,  r'(x)  = 0 is  equivalent  to 


f ( x)  = (a12-a22)x + (l-a1)(a12x2+4ck2) ^ - (l-a2)(a22x2+4ck2)iS  = 0 , (2.25) 


-10- 


where  a,  = 2(k+2)c, 


if  k > 3,  and  the  definition^  of  a 


and  a^  are  interchanged  if  k = 1.  It  will  now  be  shown  that  f as  defined  by 
(2.25)  is  a monotonically  decreasing  function,  f(0)  > 0,  f(»)  < 0,  and  hence 


the  maximum  is  the  unique  solution  of  (2.25).  By  differentiation 


Consider  k > 3,  a similar  argument  holding  for  k=l.  a < a is  equivalent 


< k/(6k+4),  true  by  Lemma  A5.  Hence  a 


By  direct  evaluation,  f(0)  > 0 and  f(<*»)  < 0,  giving  the  conclusion.  The 
maximizing  value  Xj^  can  be  conveniently  obtained  by  a very  short  FORTRAN 
program  that  does  interval  halving,  and  Table  1 was  so  obtained. 


Table  1 


Maxima  of  r(x) 


Maximum 


1000 


Suppose  k is  fixed  and  M is  the  maximum  value  of  r(x) 


obtained  as 


max 


and  the  absolute  error  for  either  approximation  is  bounded  by 
2 

(M-lKl+x  /k)f^(x)p(x,Ymjn)  and  for  any  given  x is  less  than 
2 

(1+x  /k)f,  (xKp(x,y  )-p(x,Y  . )),  and  the  absolute  value  of  the  relative 
k max  min 

error  is  bounded  by  M-l  and  for  any  given  x is  less  than  p(x,Y  )/p(x,Y  . 

r max  r mi 

Some  numerical  examples  are  given  in  Table  2.  Here  is  the  right-hand 

and  U . the  left-hand  expression  of  (2.27). 


Tabl*  2 


Upper  and  Lower  Bounds  to  t-Tail  Probabilities 


.696-01 

.696-01 

.140-01 

.140-01 

.464-02 

.464-02 

.509-01 

.512-01 

.516-02 

.517-02 

.923-03 

.924-03 

.366-01 

.369-01 

.126-02 

.126-02 

.660-04 

.662-04 

.295-01 

.298-01 

.352-03 

.353-03 

.362-05 

.365-05 

.254-01 

.257-01 

.105-03 

.105-03 

.109-06 

.110-06 

.240-01 

.243-01 

.607-04 

.611-04 

.159-07 

.159-07 

.226-01 

.231-01 

.340-04 

.3*2-04 

.138-06 

.138-08 

3.  An  Application  to  the  Calculation  of 
Bonferroni  t-Percentiles 


Let  T^,...,Tn  be  statistics  associated  with  an  experiment  such  that 
under  the  overall  null  hypothesis  HQ  each  T\  has  a t-distribution  with  k 
degrees  of  freedom,  i = l,...,n.  Then  if  it  is  desired  to  control  the 
probability  error  rate  (the  probability  of  1 or  more  false  statements  under 
H ),  one  possibility  is  to  reject  the  null-hypothesis  associated  with  T.  if 


Ti  ta/2n  * 

where  F^t^^)  = a^n»  anc*  to  acceP't  if 

lTi'  i V2n  • 1 < i 5 " • 

Here  ^a/2n  caHed  the  100  a/2n  Bonferroni  t -percentile . Similarly  we  may 
define  the  descriptive  Bonferroni  level  by 


nPHn[lTJ  > ^iJ ] ’ 


(3.1) 


where  1 1 I is  the  maximum  of  the  observed  values  of  I T, I .....  I T I,  and  the 
'm'  1 1 1 1 n1 

actual  descriptive  level,  P [Max  |t.|  > |t  |],  will  be  less  than  or  equal 

0 l<i<n  1 m 

to  the  descriptive  Bonferroni  level.  Hence  it  is  of  interest  to  find  t 

a/2n 

and  to  evaluate  probabilities  of  the  type  in  (3.1).  Recently  Bailey  (1977) 
has  given  a table  of  ta^2n  for  small  n and  n associated  with  pairwise  com- 
parisons. The  method  to  be  described  below  has  the  advantages  of  simple 
computer  implementation  (short  program,  no  tables  to  be  stored),  and  of 
working  for  any  a and  n. 

Denote  by  t the  unique  solution  in  t of 
max 


(1+t  /k)f,  (t)p(t,y  ) = a/2n  , 
k max 


(3.2) 


and  by  tm^n  the  unique  solution  in  t of 


(1+t  /k)f,  (t)p(t.y  . ) = a/2n 
K min 


Then  it  immediately  follows  from  (2.27)  that 


t . < t < t , 

min  ct/2n  max 


(3.3) 


and  hence  either  may  be  used  as  an  approximation,  or  the  digits  common  to 
both.  The  solutions  to  (3.2)  and  (3.3)  can  be  obtained  very  simply  by  a 
short  FORTRAN  program  that  does  interval  halving.  Table  3 was  generated 
this  way.  If  the  entries  in  Bailey's  Table  are  rounded  to  three  decimals, 
they  all  fall  within  the  intervals  in  Table  3.  The  maximum  difference  between 
the  upper  and  lower  bound  in  Table  3 is  .006  and  this  occurs  for  a = .05,  n = l. 
In  general,  the  smaller  a and  the  larger  n,  the  better  the  approximation. 

The  n and  a for  Table  3 were  chosen  deliberately  to  correspond  to  Bailey's 
Table. 


Table  3 

Approximations  to  Bor.ferroni  t-Pereentiles 


20  91  190  1 20  91  190 


^nin  'max  'nin  'max  'min  'max  'min  'max  'min  'max  'min  'max  'min  kmax  min  max 
3 5.641  5.842  16.326  16.327  27.131  27.131  34.698  34.698  3.182  3.184  9.465  9.465  15.816  15.817  20.253  20.253 


5 

4.032 

4.034 

7.976 

7.976 

10.962 

10.962 

12.758 

12.758 

2.570 

2.573 

5.604 

5.605 

7.817 

7.817 

9.136 

9.137 

10 

3.169 

3.172 

5.049 

5.050 

6.139 

6.140 

6.715 

6.716 

2.227 

2.231 

4.004 

4.006 

4.985 

4.988 

5.496 

5.498 

20 

2.845 

2.848 

4.146 

4.148 

4.796 

4.798 

5.116 

5.118 

2.08* 

2.090 

3.455 

3.457 

4.105 

4.107 

4.421 

4.422 

50 

2.677 

2.680 

3.723 

3.725 

4.199 

4.201 

4.423 

4.424 

2.007 

2.012 

3.184 

3.186 

3.692 

3.694 

3.928 

3.929 

100 

2.625 

2.629 

3.598 

3.600 

4.027 

4.029 

4.226 

4.227 

1.982 

1.988 

3.101 

3.104 

3.570 

3.572 

3.784 

3.785 

1000 

2.580 

2.583 

3.492 

3.494 

3.883 

3.884 

4.061 

4.062 

1.960 

1.966 

3.030 

3.033 

3.466 

3.468 

3.662 

3.664 

The  probabilities  in  (3.1)  may  be  bounded  by  (2.27)  and  numer'cal 


1 


examples  have  been  given  in  Table  2. 

The  listings  of  the  three  very  short  computer  programs  referred  to  in 
this  paper  (maximum,  bounds,  percentiles)  are  available  on  request  from  the 
author. 


> 
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Appendix:  Results  for  the  c. 


For  ease  of  reference,  all  the  necessary  results  involving  c^'s  are 
given  here  in  the  form  of  lemmas. 


Lemma  Al:  The  even  and  the  odd  c^'s  form  increasing  sequences  with  limit 

A 


1/(2tt) 


... 


Proof : That  the  limit  is  1/ ( 2tt ) is  well-known.  By  the  recursion  formula 
for  the  T function. 


C2i+3  . , i+1  ^ / 2i+l ^ 
c„_...  " (i+l/2)(2i+3)  1 • 1 - 0 


"2i+l 


Similarly 


c *5 

2i+2  _ ,i+l/2w  is..  ... 

“ (-T_)(rn:)  51 1 » 1 * 1 • 


*2i 


completing  the  proof. 


Lemma  A2:  For  k > 1,  l-4c,  > 0. 

- k 


Proof : This  is  proved  differently  in  Soms  (1977),  but  for  the  sake  of 
completeness  a shorter  proof  is  given  here.  By  Lemma  Al,  it  suffices  to 
show  that 

l/(  2tt)  < 1/4  , 

which  is  true. 


2 2 

Lemma  A3:  For  k > 3,  8c^  - 1 > 0,  and  for  k = 1,  8c^  - l < 0. 


Proof : By  direct  computation,  8c^  - 1 < 0.  From  Lemma  Al,  it  suffices 


2 — 2 

to  show  that  8c^  * 16/(n/3  ) > 1 and  8c^  = 9/8  > 1,  which  is  true. 


Lemma  A4:  For  k > 3,  < k/(2(k+2)ck2)  - 1,  and  for  k = 1 the 


inequality  is  reversed. 
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(A  I)  will  follow  from  Lemma  A5  and  (A 2)  is  verified  directly. 

2 9 

Lemma  A5:  For  k > 3,  c^  < k/(6k+4)  and  for  k = 1,  > k/(6k+4). 

Proof ; The  case  k = 1 was  shown  to  hold  in  Lemma  A4.  Since  1/2it  < l/( 6+4/16 ) , 

2 

c^.  < k/(6k+4)  for  k > 16.  It  is  verified  directly  that  it  holds  for 

3 < k < 15. 
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